nature astronomy 


Article 


https://doi.org/10.1038/s41550-022-01812-x 


Nested dust shells around the Wolf-Rayet 
binary WR140 observed withJWST 


Received: 20 August 2022 Ryan M. Lau®@' 


Accepted: 21 September 2022 


Published online: 12 October 2022 


® Check for updates 


, Matthew J. Hankins”, Yinuo Han%, loannis Argyriou ©“, 
Michael F. Corcoran®*, Jan J. Eldridge ® ’, Izumi Endo ®®, Ori D. Fox ®°, 
Macarena Garcia Marin"°", Theodore R. Gull © °"2, Olivia C. Jones ®*, 
Kenji Hamaguchi®”, Astrid Lamberts®"°, David R. Law’, Thomas Madura”, 
Sergey V. Marchenko"®””, Hideo Matsuhara”°, Anthony F. J. Moffat © 7", 
Mark R. Morris ©°, Patrick W. Morris ©”, Takashi Onaka®®”, 


Michael E. Ressler ©7°, Noel D. Richardson ©”, Christopher M. P. Russell ®°, 
Joel Sanchez-Bermudez®7°*°, Nathan Smith ®*, Anthony Soulain ® *, 
lan R. Stevens © *, Peter Tuthill**, Gerd Weigelt ©, Peredur M. Williams © *° 


and Ryodai Yamaguchi 


20,37 


Massive colliding-wind binaries that host a Wolf-Rayet (WR) star present 

a potentially important source of dust and chemical enrichment in the 
interstellar medium. However, the chemical composition and survival of 
dust formed from such systems is not well understood. The carbon-rich 
Wolf-Rayet binary WR 140 presents an ideal astrophysical laboratory for 
investigating these questions, given its well-defined orbital period and 
predictable dust-formation episodes every 7.93 years around periastron 
passage. We present observations from our Early Release Science 
programme (ERS 1349) with the James Webb Space Telescope Mid-Infrared 


Instrument (MIRI) Medium-Resolution Spectrometer and Imager that reveal 
the spectral and spatial signatures of nested circumstellar dust shells around 
WR 140. MIRI medium-resolution spectroscopy of the second dust shell and 
Imager detections of over 17 shells formed throughout approximately the 
past 130 years confirm the survival of carbonaceous dust grains from WR 140 
that are probably carriers of ‘unidentified infrared’-band features at 6.4 and 
7.7 4m. The observations indicate that dust-forming carbon-rich Wolf-Rayet 


binaries can enrich the interstellar medium with organic compounds and 
carbonaceous dust. 


Some of the first carbonaceous dust grains and organic material in 
the Universe may have been produced inthe compressed, chemically 
enriched winds of massive, evolved colliding-wind binaries hosting a 
Wolf-Rayet (WR) star of the carbon-rich (WC) subtype’ *. In the hostile 
environment around these massive and evolved stars’, dust is pro- 
duced from the compressed gas in the wind collision region where 
the stronger stellar wind from the WR star interacts with the stellar 
wind of an OB star companion® ’. Galactic colliding-wind WC binaries 


with resolvable circumstellar dust nebulae therefore provide impor- 
tant laboratories for studying this dust-formation process, where 
observations over the past few decades have demonstrated how dust 
formationis regulated by the orbit of the binary system”***. However, 
the chemical signatures of dust grains formed in colliding-wind WC 
binaries’, and the survival of these grains as they propagate into the 
interstellar medium (ISM) in their harsh radiative stellar environment, 
are not well understood. 


A full list of affiliations appears at the end of the paper. 
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Fig. 1|JWST/MIRI observations of WR 140. Left: false-colour JWST/MIRI Imager 
observations of WR 140 taken with the F770W, F1500W and F2100W filters that 
correspond to blue, green and red, respectively. Eight symmetric diffraction 
spikes are seen around the saturated core of WR 140 (at the origin) that exhibit 
bluer colours than the dust emission. Right: JWST/MIRI MRS spectral data cube 
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images of WR140 from Channel 2 at 8.01 ppm taken at two different positions 
that cover the central binary and shell 1 (Offset 1, top) and shells 1 and 2 (Offset 
2, bottom) shownina square-root stretch. The central binary and bright dust 
emission features are labelled in the MRS data cube images. 


Arguably the best example of ‘episodic’ dust formation in colliding 
winds is the WR binary WR140 (also known as HD 193793), acarbon-rich 
evolved massive star of subtype WC7 in a highly eccentric (e = 0.89), 
2,895 d (7.93 yr) mutual orbit around an early-type O5.5fc star’? ”. 
Unlike continuously dust-forming WC binaries like WR 104 (ref. °) that 
are known for their ‘pinwheel’ nebulae, WR 140 periodically produces 
dust over episodes of a few months in duration and exhibits a nebular 
morphology resembling segmented shells®. This episodic dust produc- 
tion by WR 140 occurs when the dense stellar wind of the WC7 star is 
sufficiently compressed around periastron passage by the wind of the 
Ostar®”"°. The dust condenses in the shock-compressed WC wind andis 
carried into the circumstellar medium, probably seeding the ISM with 
new carbonaceous dust. Given the relatively close distance of 1.64 kpc 
to WR 140 (refs. ”’), previous infrared (IR) imaging by ground-based 
telescopes has resolved thermal emission from (at most) two discrete, 
nested dust shells out to ~5,000 au from the central binary, and has 
provided a comparison of dust properties from one periastron pas- 
sage to the next®. However, mid-IR spectroscopy of spatially resolved 
dust emission and sensitive imaging observations detecting more 
distant dust shells are critical to investigating the chemical composi- 
tion and survival of WC dust. Such observations have been practi- 
cally impossible due to the sensitivity limitations of ground-based 
facilities and the low angular resolution of previous space-based 
observatories. 

As part of the James Webb Space Telescope (JWST) Director’s 
Discretionary Early Release Science programme (ID ERS 1349), we 
observed WR 140 withJWST to investigate the chemical composition 
and survival of dust formed by colliding winds of WC binaries. WR140 
also presented an ideal target for the Director’s Discretionary Early 
Release Science programme to demonstrate the capabilities of JWST’s 
Mid-Infrared Instrument (MIRI) to resolve and detect faint extended 
emission around a bright central point source. In this Article we discuss 


the JWST/MIRI observations that reveal the spectral and spatial signa- 
tures of dust around WR 140: over 17 nested dust shells from the past 
2130 yr of episodic dust production and the mid-IR spectroscopic 
signatures from the direct emission of WC dust. 


Results and discussion 

WR 140 was observed by JWST with the MIRI'*"° Medium-Resolution 
Spectrometer (MRS) and Imager on 2022 July 8 UT and 2022 July 27 
UT, respectively. The timing of the observations captured the binary 
system at an orbital phase of 0.7, ~5.5 yr after its last dust-formation 
episode and periastron passage on 2016 December”. These MIRI Imager 
(Fig. 1, left) and MRS (Fig. 1, right) observations provide mid-infrared 
(5-28 pm) imaging and spectroscopy of WR 140 and its circumstellar 
environment that are more sensitive than previous observations by 
over two orders of magnitude. 

The MIRI Imager observations of WR 140 at 7.7, 15 and 21 um shown 
in Fig. 1 (left) reveal thermal emission from nested and remarkably 
evenly spaced circumstellar dust shells exceeding distances of ~45” 
(-70,000 au) from the central binary. The extent of the these distant 
circumstellar shells detected around WR 140 exceeds that of all other 
known dust-forming WC systems by factors of 4 or greater”’”"*. The 
central binary and the two most recent dust shells in the Imager obser- 
vations are affected by detector saturation and/or partially overlap with 
the point spread function (PSF) of the bright core, whichis dominated 
by free-free emission from ionized stellar winds°®. However, due to the 
dispersion of incoming light, integral field unit imaging of these inner 
regions with the MRS was not affected by saturation (Fig. 1, right). 

One of the brightest repeating circumstellar dust features in the 
MIRI imaging observations is associated with the southwestern ‘Cl’ 
dust emission peak, which is labelled in the MIRI MRS spectral data 
cube image at 8.01 pm shown in Fig. 1 (right). The Cl feature was initially 
identified from ground-based mid-IR imaging of WR 140 in ref. °. The 
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Fig. 2| Views of the single-shell geometric dust model of WR 140 at the orbital 
configuration of MIRI observations. Left: dust column density in the on-sky 
projection of the geometric model, labelled with the observed dust features. 

NW, northwest; SE, southeast. Right: dust column density in the orbital-plane 
projection of the model. The dust column density shown has been normalized 

to unity by the maximum value in the each model. The cyan dot corresponds to 
the location of the central binary. The projected orbital configuration of WR 140 
in the reference frame of the WR star and the direction of motion of the O star is 
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shown in the insets of each projection. The semi-major axis of the orbit is 8.9 mas 
(14.6 au)”. The red regions in the orbital schematic around periastron passage 
correspond to true anomalies spanning —135° and +135°, where dust forms 

from the colliding winds (see ‘WR 140 geometric models’ in the Methods). Dust 
formation does not occur at a constant rate throughout periastron passage: more 
dust forms from colliding winds around the beginning and end of the orbital 
regions outlined in red”. 
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Fig. 3 | Comparison of the MIRI FISOOW image of WR 140 with the 20-shell 
geometric dust model. Left: MIRI Imager observations of WR 140 taken with the 
F1500W filter. Yellow contours show the diffraction spikes from the MIRI FISOOW 
PSF generated by WebbPSF” at levels corresponding to [10°, 10°*, 10°, 107] 
times the peak intensity of the PSF. Right: dust column density model convolved 
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with the MIRI FISOOW PSF. The dust column density shown has been normalized 
to unity by the maximum value in the model. The inset shows roughly the inner 3” 
of the dust model with the central core. The cyan dot indicates the position of the 
central binary. 


horizontal ‘Southern Bar’, which extends eastwards from C1, and the 
curved ‘Eastern Arm are also previously identified features that show 
repeating structures in the MIRI imaging observations. However, the 
Southern Bar and Eastern Arm are not detected as far out from the 
central binary as Cl. Eight prominent and symmetrical diffraction 
spikes due to the secondary struts and the hexagonal shape of JWST’s 
primary mirror are present around the saturated core of WR140 inthe 
MIRI image and exhibit bluer colours than the dust emission (Fig. 1, 
left). Interestingly, there are two additional asymmetric and ‘redder’ 
linear features that extend from the core of WR 140 to the northwest 
and southeast. These linear features are not consistent with known 
instrumental artefacts from bright point sources, which indicates that 


their origin is astrophysical and probably attributable to emission from 
circumstellar dust. 

The observed morphology of the circumstellar dust emission 
around WR 140 can be interpreted by a geometric dust shell model 
(see ‘WR 140 geometric models’ in the Methods). The models were 
generated at an orbital phase consistent with the JWST observations 
to simulate the dust column density of a single dust shell (Fig. 2) as 
well as 20 nested dust shells (Fig. 3, right). These geometric dust shell 
models utilized the well-constrained measurements of WR 140’s orbital 
parameters”’” and have been used previously to interpret the expan- 
sion of newly formed dust around periastron passage”. The single-shell 
model (Fig. 2) is shown in both on-sky and orbital-plane projections 
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Fig. 4 | Radial plots across C1 dust features from Imager observations and the 
geometric dust model. Radial profiles through the C1 dust features are shown 
for the 7.7, 15 and 21 pm (blue, green and red lines) emission and the column 
density from the dust model (black, dashed line) normalized to the peak value at 
shell 2 and offset vertically for clarity. The lower bound of the emission is clipped 
at alevel of 20. Note that the radial profiles for dust emission and column density 
do not match exactly because thermal dust emission is also regulated by dust 
temperature, which was not incorporated in the dust column density model. Grey 
vertical lines indicate intervals evenly spaced by 2.67” (4,380 au), which is the 
median separation between the flux density peaks. 


and demonstrates how column density enhancements along the line 
of sight reproduce the bright dust emission features seen in the MRS 
observations (Fig. 1, right). Specifically, the C1 feature is composed of 
dust along the projected southwestern edge of the shell. Similarly, the 
Eastern Arm and Southern Bar arise due to enhanced column densities 
along the projected eastern and southern edges of the shell, respec- 
tively. Note that when WR 140 forms dust, it does not doso at aconstant 
rate’. Enhanced dust formation occurs from colliding winds around the 
beginning and end of periastron passage, which explains the relative 
absence of dust to the northwest of WR 140 (Figs. 1 and 2, see ‘WR 140 
geometric models’ inthe Methods). 

The 20-shell geometric model presents a striking match to the 
overall morphology of WR 140’s nested dust shells revealed by the 
MIRI Imager (Fig. 3). The asymmetric and linear northwest and south- 
east features are present in the model at an angle and radial extent 
consistent with the linear features from the Imager observations. A 
comparison to a model MIRI FISOOW PSF generated by WebbPSF”° 
also demonstrates that the linear features do not align with the eight 
diffraction spikes (Fig. 3, left). Given the orientation of the projected 
northwestern and southeastern edges of the single dust shell (Fig. 2, 
left), the linear northwest and southeast dust features are probably 
associated with the projection of those edges from repeating dust 
shells. The model also closely reproduces the repeating bright dust 
features, as well as the regular spacing between the dust shells. A radial 
plot of the dust column density from the model and the observed dust 
emission in the direction of C1 from shell 2 and beyond are shown in 
Fig. 4 and demonstrate the consistency in the shell spacing. 

The regular spacing of the observed emission from 17 shells and 
the agreement with the geometric model illustrate an apparently 
freely expanding trajectory of the dust shells and WR140’s consistent 


Table 1 | Projected distances of C1 features from WR 140 


Shell Tel") Arg") T¢,(au) Ar,,(au) 
1 1.63 - 2,670 - 
2 4.32 2.69 7,090 4,410 
3 6.98 2.66 11,440 4,360 
4 9.62 2.64 15,780 4,330 
5 12.33 271 20,220 4,440 
6 15.01 2.68 24,610 4,390 
7 17.66 2.65 28,960 4,350 
8 20.24 2.58 33,200 4,240 
9 23.08 2.84 37,850 4,650 
10 25.76 2.68 42,250 4,400 
1 28.41 2.65 46,590 4,340 
a2, 31.14 213 51,070 4,480 
3 33.78 2.64 55,400 4,330 
14 36.56 2.78 59,960 4,560 
5 39.18 2.62 64,250 4,300 
16 41.98 2.80 68,850 4,590 
7 44.60 2.61 73/40 4,290 


Io, is the projected separation distance between C1 features and WR 140, which was 
calculated from the F1500W image for shells 2-17 and the MRS observation for shell 1. Arg, is 
the separation distance with respect to the adjacent inner dust shell adopting a distance of 
1.64 kpc (refs. ”"°). 


episodic dust-production history over approximately the past 130 yr. 
The median separation of the first 17 dust shells along Clis 2.67 + 0.07” 
(Table 1), which corresponds to 4,380 + 120 au assuming a distance 
of 1.64 kpc to WR 140 (refs. *"*). The separation distance matches 
the expected separation given the 326 mas yr‘ (2,540 kms‘) proper 
motion of Cl (ref. °) over WR140’s 7.93 yr orbital period. The projected 
expansion velocity of the dust shells calculated from the median shell 
separations and the orbital period is ~-2,600 kms“, which is consist- 
ent with the previously measured C1 proper motion. As the geomet- 
ric models indicate that the expansion direction of the Cl feature in 
three-dimensional space is nearly perpendicular to our line of sight, 
the projected velocity of C1 should be consistent with the true dust 
expansion velocity”. The dust expansion velocity of -2,600 kms7is 
therefore slightly less than the terminal velocity of the WC-star wind 
(v.,= 2,860 km s* ref. 74). 

The apparently constant expansion velocity of the dust shells is 
notable given that dust around WR 140 should experience acceleration 
and deceleration due to radiation pressure and drag forces, respec- 
tively®’. Although beyond the scope of this work, the dust dynamics 
and the balance between radiation pressure and drag forces in the 
surrounding environment of WR 140 warrant further investigation. 
Given that the distant shells approach the signal-to-noise detection 
limits and retain a relatively constant separation distance (Fig. 4 and 
Table 1), fainter and cooler dust shells probably exist beyond those 
detected by the MIRI Imager observations. The detection of distant 
shells from dust production ~-130 yr ago and the expected presence 
of more distant shells indicates that the dust formed by WR 140 can 
survive its harsh circumstellar environment and probably enriches the 
surrounding ISM. Follow-up observations with JWST of the local ISM 
and the surrounding -10 pc bubble that may have been carved out by 
WR140 (ref. ”) could provide valuable insight into the ISM enrichment. 

Understanding the chemical properties and spectral signatures 
of dust formed by WC binaries like WR 140 is important given their 
potential role as dust sources in the ISM (for example, refs. 1777). 


Nature Astronomy 


Article 


https://doi.org/10.1038/s41550-022-01812-x 


20,000 + 
or 1,400.00 
18,000 4 
1,200.00 
16,000 4 
1,000.00 
ie, POS S 14000 4 | 
Ly = 3 
3 gooo0 = = [Sv] 10.5 um | 
2 a. o  % 12,0004 
= & f=) 
fo) etl 7 \ 60000 - & Artefact 
j ! 10,000 4 
Se yf 
aaa 400.00 
8,000 + 
Re 200.00 UIR-band 
a 6,000 — features 
(0) 
T T T 4,000 T T T T T T T T i 
0 -2.0 -4.0 6 ri 8 9 10 11 12 13 14 
Offset (") Wavelength (um) 
12 
1.0 4 
0.8 + 
x 
P| 
= 
B os 
NN 
o 
E 
S o4- 
02 4 
ie) eS 
T T T T T T 
6.0 6.5 7.0 75 8.0 8.5 9.0 


Wavelength (um) 


Fig. 5 | MIRI MRS spectroscopy of the C1 feature in shell 2. Top left: MRS 
spectral data cube image of WR 140 at 8.01 ppm (Channel 2) taken at the Offset 

2’ position overlaid with the apertures used to extract the spectra from the C1 
feature in shell 2 (solid red circle) and measure the background emission (dashed 
red circle). Top right: extinction-corrected and background-subtracted 6-14 pm 
spectra of the Cl feature in shell 2 (red line) overlaid with the continuum fit froma 
cubic Chebyshev polynomial function (black line). The absorption ‘feature’ near 


12 pmis an instrumental artefact, and the emission feature at 10.5 mis probably 
emission from the [S Iv] 10.5 jum fine-structure line. Bottom: 6-9 pm continuum- 
subtracted and normalized spectrum of the Cl feature from shell 2 exhibiting 
prominent 6.4 and 7.7 1m UIR-band features. Gaussian-fitted peaks of these 
features are located at 6.38 um and 7.75 ppm and are indicated by dashed vertical 
grey lines. Spectra were smoothed by convolving with a box filter kernel witha 
width of 21 spectral elements. 


We can now investigate this using spatially and spectrally resolved 
observations of WR 140’s circumstellar dust with the MIRI MRS. The 
6-9 um wavelength interval includes a wealth of features that trace 
species of interstellar molecules and dust (see ref. **) such as polycy- 
clic aromatic hydrocarbons, which are believed to be carriers of the 
unidentified infrared (UIR) bands most commonly observed at 3.3, 
6.2, 7.7, 11.2 and 12.7 pum (refs. **”’). Polycyclic aromatic hydrocarbons 
are known to be highly stable due to the honeycomb-like structure 
of carbon atoms fused in aromatic rings”””*. The survival of WR 140’s 
circumstellar dust population in the harsh radiation field around the 
central binary would therefore be consistent with a composition of 
C-rich aromatic compounds. The presence of UIR bands that peak at 
6.4 and 7.9-8.5 um has been reported from other dust-forming WC 
stars’*”°; however, previous spectroscopic observations were unable 
to spatially resolve circumstellar dust emission that blended witha 
bright central core. Mid-IR integral field unit spectroscopy of WR 140 
with MIRI MRS can therefore be used to confirm the association of such 
UIR features with the WC dust. 


Figure 5 presents the spectral data cube image of WR 140 at 
8.01 pm, the extinction-corrected spectrum of C1 from shell 2 along 
with the dust-continuum fit (see ‘Dust-continuum subtraction’ in the 
Methods), and the 6-9 pm dust-continuum-subtracted and normal- 
ized spectrum from the shell 2 C1 feature. We identified the presence 
of two prominent UIR bands at 6.4 and 7.7 pm, which resemble the 
features from previous unresolved spectroscopic observations of 
other dust-forming WC stars. The peak positions of these features 
appear at 6.38 pm and 7.75 um according to a simple Gaussian fit to the 
continuum-subtracted spectrum (Fig. 5, bottom). The 7.7 pm feature is 
typically attributed to aromatic C-C stretching modes in polycyclic aro- 
matic hydrocarbons”. The 6.4 tim feature, whichis primarily detected 
in H-poor environments from episodic dust-forming sources such as 
novae’’, R Coronae Borealis stars” and other WC stars’, is thought to 
originate from large carbonaceous molecules or small, H-free carbona- 
ceous grains”’*’. The detection of the 6.4 um band from circumstellar 
dust around WR 140 is therefore consistent with its expected H-poor 
environment dominated by the C-rich WC-star wind. The common 
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Table 2|WR140 MIRI observation summary 


MIRI MRS MIRI Imager 


Coordinates (J2000) 20h 20min 27.985s, 2h20min27.976s, 


+43°51'15.90" +43°51'16.28" 
20h20 min27.781s, 
+43°51'13.05" 
Observation date (UT; 2022 July 8; MJD 59768 2022 July 27; MJD 
MJD) 59787 
WR 140 phase (¢) 0.70 0.71 
Channels/filters Channels 1- 4 F770W, F1500W, 
F2100W 
Wavelength (um) Ch 1: 4.9-7.65 7.7, 15.0, 21.0 
Ch 2: 7.51-11.71 
Ch 3: 11.55-18.03 
Ch 4: 17.71-28.1 
Spectral resolving power 1,500-3,500 0.25, 0.48, 0.67 
(A/AA) 
Spaxel/pixel scale (") 0.13-0.35 O11 
Field of view 3.2"*2.7"-6.6"*7.6" 112.6"*73.5" 


The two sets of MRS coordinates correspond to the pointings of the Offset 1 and Offset 2 
observations (Fig. 1, right), and the Imager coordinate indicates the pointing of the Imager 
observation. @ is WR 140’s orbital phase at the time corresponding to the MRS and Imager 
observations based on the orbit described in ref.'*. The integral field unit channels and 
wavelength ranges used in the MRS observations are provided, as well as the Imager filters 
and their corresponding central wavelengths. 


8.6 pum UIR feature, which is attributed to C-H bending modes, is nota- 
bly absent from the spectrum and is also consistent with an H-poor 
environment. The direct detection of UIR bands from shell 2 indicates 
that dust formed by WR140 is carbonaceous and richin aromatics. The 
stability of such aromatic compounds””* supports the interpretation 
of WC dust survival out to the distant shells, as revealed by the Imager 
observations (Figs. land 4). 

Due to their prevalence throughout the Universe and their influ- 
ence on the thermal regulation and chemistry of the ISM (see ref. * 
and ref. therein), carriers of these UIR-band features are important 
components of the ISM. The dominant source of UIR-band carriers 
has widely been thought to be evolved, solar-mass asymptotic giant 
branch stars”***, However, the timescale for massive stars to evolve 
to the WR phase (millions of years) is much shorter than the time- 
scale for lower-mass stars to evolve into asymptotic giant branch stars 
(2100 Myr). On the basis of a dust-production analysis with Binary 
Population and Spectral Synthesis® models in ref. ’, the population of 
WC binaries could produce as much dust as the population of asymp- 
totic giant branch stars, depending on star-formation histories and 
metallicities. The detection of the Cl feature out to distant shells (for 
example shell 17; Fig. 4) also demonstrates that the carbonaceous WC 
dust persists in the luminous and hard radiation field of the central 
binary system for at least 130 yr after the initial dust-formation event 
and can probably propagate to the surrounding ISM. We therefore 
argue that dust-forming WC binaries like WR 140 could be considered as 
apossible early and potentially dominant source of organic compounds 
and carbonaceous dust in the ISM of our Galaxy and galaxies beyond. 


Methods 

JWST/MIRI observations 

MIRI Imager. WR 140 (J2000 RA 20 h20 min 27.9 s, dec. +43° 51’ 16.29”; 
ref. °°) was observed with the JWST/MIRI Imager on 2022 July 27 UT with 
the F770W, F1500W and F2100W filters inthe FULL subarray mode. The 
Imager observations in each filter were performed using three sets of 
four-point dithers at starting set 5 in the positive direction with the 


default pattern size that was optimized for an extended source. The 
FASTRI1 readout pattern with 43 groups, 1 integration and 1 exposure 
per dither was used for each filter observation. The total exposure time 
in each filter was 0.40 h. Asummary of the MIRI Imager observations 
of WR 140 is provided in Table 2. 

The MIRI Imager data were processed with the JWST calibration 
pipeline version 1.6.2 and STCal version 1.0.0. After running the detec- 
tor level pipeline, the rate file was recreated from the rateints file to 
ensure that the correct rate value was reported in regions that saturated 
after the third group. We note that a dedicated nonlinearity correction 
for the F2100W data was not available when we reduced this dataset. As 
the MIRI detectors exhibit a wavelength dependency in their nonlinear 
behaviour, dedicated corrections were needed at wavelengths redder 
than 18 pm. We therefore expected an additional error contribution of 
a few per cent at 21 um. 


MIRI medium-resolution spectroscopy. WR 140 was observed with 
the MIRI MRS” on 2022 July 8 UT. MIRI MRS observations provided 
spatially resolved spectroscopy over the four channels (Ch 1-Ch 4) 
that covered the 4.9-28.1 pm interval with a spectral resolving power of 
R=1,500-3,500 (ref. *’”). Two sets of MRS observations that targeted the 
two most recent dust shells formed by WR 1405.5 and 13.4 yr ago were 
performed (Fig. 1, right). The coordinates of the MRS observations and 
the other observing details are provided in Table 2. The two MRS obser- 
vations were linked in a non-interruptible sequence due to the high 
proper motion of the circumstellar dust shells (-300 mas yr“; ref. °). 

The first MRS observation utilized the bright core of WR 140 for 
target acquisition, which used the neutral density filter for 11.1sinthe 
FAST readout pattern. Target acquisition was not used for the second 
observation in the non-interruptible sequence given the small offset 
(-4”) required from the initial position. Each MRS observation was per- 
formed using the four-point dither optimized for an extended source 
and oriented in the negative direction. 60 groups and 5 integrations 
were used for each dither position, corresponding to 0.94 hof exposure 
time in each grating position (SHORT, MEDIUM, LONG). Simultaneous 
imaging with the MIRI Imager was also utilized in the MRS observations 
to improve the astrometric accuracy of the MRS data. 

Calibration level 1 MIRI MRS data of WR 140 were obtained from 
the Mikulski Archive for Space Telescopes (MAST). These level 1 data 
were then processed throughJWST calibration pipeline version 1.6.0, 
where ‘output_type = band’ was set in the ‘cube_build’ step to output 
12 data cubes corresponding to the three grating positions of each of 
the four channels. 


Data analysis 
MIRI Imager background subtraction, radial dust emission profile 
and C1 positions. The median background emission in relatively empty 
patches of sky beyond the circumstellar shells in the F770W, FISOOW 
and F2100W MIRI Imager observation exhibited fluxes of 10.4, 41.4 and 
210.2 MJy sr“, respectively. A background subtraction using the median 
values was therefore applied to the three MIRI Imager observations. 
A radial 3-pixel-width line-cut plot originating at the position of 
WR 140 and oriented 53° south of west was used to measure the emis- 
sion and positions of the C1 features from the MIRI Imager observa- 
tions (Fig. 4; Table 1). A Gaussian profile was used to measure the peak 
positions of each of shells 2-17 from the radial dust emission profile of 
the FISOOW image. An identical line cut was used to measure the dust 
column density from the 20-shell geometric dust model. 


MIRI MRS spectroscopy of C1. Spectra were extracted from the 
MIRI MRS ‘Offset 1’ observations of Clin shell 1(20 h 20 min 27.838 s, 
+43° 51’ 14.69”) using a circular aperture with a radius r= 0.7”. Sub- 
traction of the background emission and contamination from the dif- 
fraction pattern of the bright core was performed by using a circular 
aperture in a relatively dust-free position northwest of WR 140 (20h 
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20 min 27.932 s, +43° 51’ 17.75”) at the same separation distance from 
the Cl aperture. However, contamination from emission lines from 
the bright stellar core seemed to persist in the shell 1 Cl spectra even 
after background subtraction. Follow-up work from the WR Dust- 
ERS programme (ERS 1349) will be conducted to implement a tech- 
nique to perform PSF subtraction on the MIRI MRS observations ina 
future study. 

Acircular aperture of r= 0.7” was also used to extract the spectra 
from the MIRI MRS ‘Offset 2’ observations of Clin shell 2 (20 h20 min 
27.657 s, +43° 51’ 12.970”). Background emission for the C1 region in 
shell 2 was measured and subtracted using a circular background 
aperture with an identical radius centred on arelatively empty areain 
the field of view to the southwest of shell 2 at the following coordinates: 
20 h20 min 27.575 s, +43° 51’ 11.852”. The 6.4 and 7.7 um features were 
present in background-subtracted spectra even when using different 
background positions in other relatively empty regions of the sky. 

To mitigate high-frequency noise from artefacts suchas ‘fringing’, 
the MRS spectra were smoothed by convolving with a box filter kernel 
that had a width of 21 spectral elements (AA = 0.04 pm). The width of 
the UIR features (20.2 um) was notably larger than the width of the 
kernel used for smoothing. 


Spectral stitching between MRS channels. To combine the spectra 
extracted from the three subbands from channels 1-4, the median flux 
density in the overlapping wavelengths between the 12 subbands was 
used to combine the extracted Cl spectra. The spectra were normalized 
to the Channel 2 SHORT subband (A = 7.51-8.76 pm), which overlaps 
with the 6-9 pm wavelength range studied in this work. 


Dust-continuum subtraction. The dust-continuum emission was fitted 
and subtracted from the ~6—9 pm wavelength range of the extracted 
spectra of the Cl region from shell 2 using a cubic Chebyshev polyno- 
mial function from the ‘Spectrum1D’ and ‘fit_continuum’ functions 
in the Python package specutils**. Wavelengths between 6.8-7.2 1m 
and 9-15 pm were used to fit the C1 continuum in shell 2. lonized gas 
probably attributable to the [S Iv] 10.5 um fine-structure line was also 
detected in the C1 spectra and will be investigated in future analysis. 
The apparent absorption ‘feature’ near 12 pm was an instrumental 
artefact. The shorter wavelengths (A < 6 pm) were omitted for the shell 
2 dust-continuum fit to avoid enhanced emission probably arising from 
ahotter dust component (R.M.L. et al., manuscriptin preparation). The 
Cl spectra from shell 2 and the dust-continuum fit are shown in Fig. 5 
(top right). The ~6-14 pm continuum-subtracted spectra were then nor- 
malized to the peak flux density of the 7.7 1m feature (Fig. 5, bottom). 


Extinction correction. Extracted MIRI MRS spectra from the circum- 
stellar dust around WR 140 were dereddened using the extinction 
correction factors derived from the Galactic interstellar extinction 
curve from Gordon et al.*’, where the optical total-to-selective 
extinction ratio (Ry) is 3.17 and the extinction at the K,-band is 
Ax, = 0.24 (ref. ). 


WR 140 geometric models 

The circumstellar dust structures around WR 140 shown in Fig. 3 (left) 
were simulated using the geometric model presented in ref. ’. The 
model assumes dust to be produced ona thin, conical surface originat- 
ing from the wind-shock interface and has been applied in previous 
studies of dust-forming colliding-wind binaries’”*°"*. The geometry 
of the dust structure depends in part on the orbit of the binary, for 
which we adopted the well-constrained orbital parameters from ref. ”. 
The half opening angle of the conical dust surface was assumed to be 
40° (refs. ®'*“"), and the dust expansion velocity was assumed to be 
2,600 km s“, which was consistent with the velocity measured from the 
separation distance between the dust shells (Table 1). The only parame- 
ter of the geometric model that deviated from the original model fitted 


to ground-based observations of a single shell from ref. ’’ was the dust 
expansion speed, which was slightly larger than the expansion speed 
of 2,450 kms“ fitted to data at an earlier orbital phase (0.592). This 
apparent discrepancy may reflect acceleration during the expansion of 
the dust plume found by ref. '°. The model was then convolved with the 
F1500W MIRI Imager PSF from the WebbPSF Python package”’. Note 
that because only the dust morphology was generated by the model, 
the central binary itself did not appear. 

Observations of WR 140 indicate that dust does not form evenly 
over the orbit®”’. By fitting to the structures in one shell of dust 
observed by ground-based AO-assisted IR imaging, ref. ’° found that the 
structures were reproduced if dust production from the colliding winds 
began at a true anomaly of -135° and ceased at +135° (Fig. 2), between 
which dust production was suppressed. No dust was assumed to be pro- 
duced along other parts of the orbit. Furthermore, the structures were 
better reproduced when it was assumed that the dust density peaked 
onthe trailing edge of the conical shock front, becoming progressively 
lower when moving towards the leading edge. Further details of the 
geometric model of WR140 are provided in ref. ”. 


Data availability 

Data used in this study were obtained under JWST Director’s Discretion- 
ary Early Release Science programme ID ERS 1349 and have no exclusive 
access period. Datacan be obtained from the Mikulski Archive for Space 
Telescopes (MAST; https://archive.stsci.edu/missions-and-data/jwst). 


Code availability 

This research made use of Astropy, a community-developed core 
Python package for Astronomy”. This research also made use of Jdaviz, 
which is a package of astronomical data analysis visualization tools 
based onthe Jupyter platform*’. 
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